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Using the well-known "displace, cut and reflect" method used to generate disks from given solu- 
tions of Einstein field equations, we construct static disks made of perfect fluid based on vacuum 
Schwarzschild's solution in isotropic coordinates. The same method is applied to different exact 
solutions to the Einstein's equations that represent static spheres of perfect fluids. We construct 
several models of disks with axially symmetric perfect fluid halos. 

All disks have some common features: surface energy density and pressures decrease monotonically 
and rapidly with radius. As the "cut" parameter a decreases, the disks become more relativistic, 
with surface energy density and pressure more concentrated near the center. Also regions of unstable 
circular orbits are more likely to appear for high relativistic disks. Parameters can be chosen so 
that the sound velocity in the fluid and the tangential velocity of test particles in circular motion 
are less then the velocity of light. This tangential velocity first increases with radius and reaches a 
maximum. 
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I. INTRODUCTION 



Axially symmetric solutions of Einstein's field equations corresponding to disklike configurations of matter are 
of great astrophysical interest, since they can be used as models of galaxies or accretion disks. These solutions can be 
static or stationary and with or without radial pressure. Solutions for static disks without radial pressure were first 
studied by Bonnor and Sackfield and Morgan and Morgan |2j, and with radial pressure by Morgan and Morgan 
[3 . Disks with radial tension have been considered in Q| , and models of disks with electric fields [I| , magnetic fields 
[3 , and both magnetic and electric fields have been introduced recently Q ■ Solutions for self-similar static disks were 
analyzed by Lynden-Bell and Pineault [a, and Lemos M\. T he superposition of static disks with black holes were 
considered by Lemos and Letelier [Tol llll IT ^ |. and Klein|i1| . Also Bicak, Lynden-Bell and Katz Q studied static 
disks as sources of known vacuum spacetimes and Bicak, Lynden-Bell and Pichon |l5j found an infinity number of 
new static solutions. For a recent survey on relativistic gravitating disks, see (lt^ . 

The principal method to generate the above mentioned solution is the "displace, cut and reflect" method. One 
of the main problem of the solutions generated by using this simple method is that usually the matter content of 
the disk is anisotropic i.e., the radial pressure is different from the azimuthal pressure. In most of the solutions the 
radial pressure is null. This made these solutions rather unphysical. Even though, one can argue that when no radial 
pressure is present stability can be achieved if we have two circular streams of particles moving in opposite directions 
(counter rotating hypothesis, see for instance 

In this article we apply the "displace, cut and reflect" method on spherically symmetric solutions of Einstein's 
field equations in isotropic coordinates to generate static disks made of a perfect fluid, i.e., with radial pressure equal 
to tangential pressure and also disks of perfect fluid surrounded by an halo made of perfect fluid matter. 

The article is divided as follows. Section II gives an overview of the "displace, cut and reflect" method. Also 
we present the basic equations used to calculate the main physical variables of the disks. In section III we apply 
the formalism to obtain the simplest model of disk, that is based on Schwarzschild's vacuum solution in isotropic 
coordinate. The generated class of disks is made of a perfect fluid with well behaved density and pressure. Section IV 
presents some models of disks with halos obtained from different known exact solutions of Einstein's field equations 
for static spheres of perfect fluid in isotropic coordinates. In section V we give some examples of disks with halo 
generated from spheres composed of fluid layers. Section VI is devoted to discussion of the results. 
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FIG. 1: An illustration of the "displace, cut and reflect" method for the generation of disks. In A the spacetime with a 
singularity is displaced and cut by a plane (dotted line), in B the part with singularities is disregarded and the upper part is 
reflected on the plane. 

II. EINSTEIN EQUATIONS AND DISKS 

For a static spherically symmetric spacetime the general line element in isotropic spherical coordinates can be 
cast as, 

ds 2 = e v ^dt 2 - e A « [dr 2 + r 2 (d9 2 + sin 2 Od^ 2 )} . (1) 
In cylindrical coordinates (t,R,z,ip) the line element takes the form, 

ds 2 = e v{R,z) dt 2 _ e X{R,z) ( di? 2 + d2 2 + R 2 d(p 2^ (2) 

The metric of the disk will be constructed using the well known "displace, cut and reflect" method that was 
used by Kuzmin |17| in Newtonian gravity and later in general relativity by many authors p|-|l6|]. The material 
content of the disk will be described by functions that are distributions with support on the disk. The method 
can be divided in the following steps that are illustrated in Fig. ^ First, in a space wherein we have a compact 
source of gravitational field, we choose a surface (in our case, the plane z = 0) that divides the space in two pieces: 
one with no singularities or sources and the other with the sources. Then we disregard the part of the space with 
singularities and use the surface to make an inversion of the nonsingular part of the space. This results in a space 
with a singularity that is a delta function with support on z = 0. This procedure is mathematically equivalent to 
make the transformation z — » \z\ + a, with a constant. In the Einstein tensor we have first and second derivatives 
of z. Since d z \z\ = 2d(z) — 1 and d zz \z\ = 2S(z), where $(z) and S(z) are, respectively, the Heaviside function and 
the Dirac distribution. Therefore the Einstein field equations will separate in two different pieces [l^: one valid for 
z 7^ (the usual Einstein's equations), and other involving distributions with an associated energy- momentum tensor, 
Tat = QabS(z), with support on z = 0. For the metric l|2|l. the non-zero components of Q a b are 

Q\ = ^[-b zz +9 zz K + bl+b%)], (3) 
QR l = Q% = ^[-b zz +9 zz {b\ + b* + bl)], (4) 
where b a b denote the jump of the first derivatives of the metric tensor on the plane z = 0, 

bab = 9ab,z\z=0+ - 9ab,z\z=0- > (5) 

and the other quantities are evaluated at z = + . The "true" surface energy- momentum tensor of the disk can be 
written as S a b — y/—g Z zQab, thus the surface energy density a and the radial and azimuthal pressures or tensions (P) 
read: 

<r = V=ihMt, P = -V=gTzQ% = -V z: 9TzQ l £- (6) 

Note that when the same procedure is applied to an axially symmetric spacetime in Weyl coordinates we have = 0, 
i.e., we have no radial pressure or tension. 

This procedure in principle can be applied to any spacetime solution of the Einstein equations with or without 
source (stress tensor). The application to a static sphere of perfect fluid is schematized in Fig. [21 The sphere is 
displaced and cut by a distance a less then its radius r&. The part of the space that contains the center of the sphere 



FIG. 2: An illustration of the "displace, cut and reflect" method for the generation of disks with halos. In A the sphere of 
perfect fluid is displaced and cut by a plane (dotted line), in B the lower part is disregarded and the upper part is reflected on 
the plane. 




FIG. 3: An illustration of the "displace, cut and reflect" method for the generation of disks with various layers of halos. In A 
the sphere with different layers of fluid is displaced and cut with a plane (dotted line), in B the field is reflected on the plane. 

is disregarded. After the inversion of the remaining space, we end up with a disk surrounded by a cap of perfect fluid. 
The properties of the inner part of the disk will depend on the internal fluid solution, but if the internal spherical 
fluid solution is joined to the standard external Schwarzschild solution, the physical properties of the outer part of 
the disk will be those originated from Schwarzschild's vacuum solution. 

In isotropic coordinates the matching at the boundary of the fluid sphere leads to four continuity conditions 
metric functions e A and e v together with their first derivatives with respect to the radial coordinate should be 
continuous across the boundary. In addition, to have a compact body the pressure at the surface of the material 
sphere has to drop to zero. Also to have a meaningfully solution the velocity of sound, V 2 = * should be restricted 
to the interval < V < 1. 

The Einstein equations for a static spherically symmetric space time in isotropic coordinates for a perfect fluid 
source give us that density p and pressure p are related to the metric functions by 



where primes indicate differentiation with respect to r. 

Also static spheres composed of various layers of fluid can be used to generate disks with halos of fluid layers 
(see figure OJl . The disk will then be composed of different axial symmetric "pieces" glued together. The matching 
conditions at the boundary of adjacent spherical fluid layers in isotropic coordinates involves four continuity conditions: 
the two metric functions e A and e v , the first derivative of A with respect to the radial coordinate, and the pressure 
should be continuous across the boundary. At the most external boundary, the metric functions e A and e v , and their 
first derivatives with respect to the radial coordinate should be continuous across the boundary; also the pressure 
there should go to zero. 



P 




(7) 



P = e - i( AT n AV + -(YW) 



(8) 
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III. THE SIMPLEST DISK 

We first apply the "displace, cut and reflect" method to generate disks discussed in the previous section and 
depicted in Fig. ^ to the Schwarzschild metric in isotropic coordinates (t,r,0,ip), 

-) 2 ^ {, , m Y 



d.s 2 = V ^ 2 dt 2 -(l + ^-) [dr 2 + r 2 (d0 2 +sin 2 0d^ 2 )] . (9) 

Expressing solution © in cylindrical coordinates, and using Eq. |J3j ©j we obtain a disk with surface energy density 
a and radial and azimuthal pressures (or tensions) P given by 



4ma 



ir(m + 2VR 2 +a 2 ) 3 ' 



(10) 



2m 2 a , 

p = ; ; 11 

7r(m + 2Vi? 2 + a 2 ) 3 (m - 2y/R 2 + a 2 ) 

The total mass of the disk can be calculated with the help of Eq. 11(J|) : 

r°° r 2lT m 

M= / o^Jg RR g w dRdip = — (m + 4a). (12) 

Eq. (|1U|I shows that the disk's surface energy density is always positive (weak energy condition). Positive 
values (pressure) for the stresses in azimuthal and radial directions are obtained if m < 2 V R 2 + a 2 . The strong 
energy condition, a + P vv + P RR = a + 2P > is then satisfied. These properties characterize a fluid made of matter 
with the usual gravitational attractive property. This is not a trivial property of these disks since it is known that the 
"displace, cut and reflect" method sometimes gives disks made of exotic matter like cosmic strings, see for instance 

Another useful parameter is the velocity of sound propagation V, defined a,sV 2 — which can be calculated 
using Eq. (JTUJ and Eq. ijTTjl: 



r2 m(4v R 2 + a 2 — m) 



v = v v ; ==?-. (13) 

3(m-2Vi? 2 + a 2 ) 2 



The condition V 2 < 1 (no tachyonic matter) imposes the inequalities m < V R 2 + a 2 or m > 
pressure condition and the speed of sound less then the speed of light condition are to be simultaneously satisfied, 
then m < \/R 2 + a 2 . This inequality will be valid in all the disk if m < a. 

With the presence of radial pressure one does not need the assumption of streams of rotating and counter 
rotating matter usually used to explain the stability of static disk models. However, a tangential velocity (rotation 
profile) can be calculated by assuming a test particle moves in a circular geodesic on the disk. We tacitly assume that 
this particle only interact gravitationally with the fluid. This assumption is valid for the case of a particle moving in 
a very diluted gas like the gas made of stars that models a galaxy disk. 

The geodesic equation for the R coordinate obtained from metric is 

e*R + i(e" W 2 - \y). R (R 2 + i 2 )-\ (e X R 2 ), R ^ = 0. (14) 



For circular motion on the z — plane, R = R = and z = 0, then Eq. (|14(l reduces to 



^ in. 



R 



i 2 {e x R 2 ),R 

The tangential velocity measured by an observer at infinity is then 

2 _ 9 VV ( d(f\ 2 _ „ 2 e A (e iy ) ! 



From the metric on the disk, 



(15) 



2<JR 2 +a 2 



aU<i e " = < :l + 2VR 2 1 ' M?i 
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FIG. 4: (a) The surface energy density a, (b) pressures P, (c) sound velocity V and (d) tangential velocity v c (rotation curve 
or rotation profile) with m = 0.5 and a = 0.6, 0.8, 1.0, and 1.2 as function of R = R/m. We use geometric units G = c = 1. 



we find that Eq. lfTp|l can be cast as, 

v? = . s ^ • (18) 

l -wm^) [(^ 2 + « 2 ) 3/2 + f(« 2 -^ 2 )] 

For R >> a, Eq. (|18fl goes as v c = (m/i?) 1 / 2 , the Newtonian circular velocity. 

To determine the stability of circular orbits on the disk's plane, we use an extension of Rayleigh [23, H]] criteria 
of stability of a fluid at rest in a gravitational field 

^>0, (19) 

where h is the specific angular momentum of a particle on the disk's plane: 

dip dp dt 

h = - 9 ^Ts = - 9 ^Ts- (20) 

Using Eq. (|15|) and the relation 

/dt Y p 2„A/^^ Xi 

one obtains the following expression for /i: 



1 = e "U ~ Re ■ (21 ) 



* ~ ^ V e"(i? 2 e A ),l -ik^e"),* ' (22) 



6 



0.035 ■ 



0.030 - 

1 0.025 

0.020 - 

0.015 - 

0.010 - 

0.005 - 

D — 




••• ::: "»»;?h«i» 



iiuimm,,, 



°f"~''- •f'.llllli 



m = 0.2 
m == 0.4 
m == 0.6 
m = 0.8 



FIG. 5: (a) The surface energy density a, (b) pressures P and (c) sound velocity V with a — 1.0 and m = 0.2, 0.4, 0.6, and 0.1 
as function of 7? = R/a. 

For the functions {I7|), Eq. reads 

2^i? 2 (l + 5;fe )V + a 2 ) 1/4 



^(i? 2 + a 2 ) 2 - 8mR 2 VR 2 + a 2 + m 2 (R 2 - a 2 ) 



(23) 



The stability criterion is always satisfied for — > 1.016. 

In figurcg](a)-(d) we show, respectively, the surface energy density, pressures, the sound velocity and curves of 
the tangential velocity (rotation curves) Eq. i|18|) with m = 0.5 and a = 0.6, 0.8, 1.0, and 1.2 as functions of R = R/m. 
Figure[S](a)-(c) display, respectively, the surface energy density, pressures and sound velocity with parameters a = 1.0 
and m = 0.2, 0.4, 0.6, 0.8 as as functions of R = R/m. We see that the first three quantities decrease monotonically 
with the radius of the disk, as can be checked from Eq. H10J1 , Ulljl and (jlHJI . Energy density decreases rapidly enough 
to, in principle, define a cut off radius and consider the disk as finite. 

IV. DISKS WITH HALOS 

Now we study some disks with halos constructed from several exact solutions of the Einstein equations for 
static spheres of perfect fluid. A survey of these class of solutions is presented in [23] • 



Buchdahl's Solution 



The first situation that we shall study is similar to the one depicted in Fig. wherein we start with a sphere of 
perfect fluid. This case will not be exactly the same as the one presented in the mentioned figure because the sphere 
has no boundary. Hence the generated disk will be completely immersed in the fluid. An example of exact solution 
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of the Einstein equations that represent a fluid sphere with no boundary is the the Buchdahl solution that may be 
regarded as a reasonably close analogue to the classical Lane-Emden index 5 polytrope |23j . The metric functions for 
this solution are: 



1 



Vl+kr 2 



e A = 1 



A 



\A + kr 2 



(24) 



where A and k are constants. Far from the origin the solution goes over into the external Schwarzschild metric, when 



The density, pressure and sound velocity are given by: 

3Ak 



kA 2 



P = 



2n(A + VI + kr 2 ) 5 ' P 2tt(-A + y/l + kr 2 ){A + y/l + kr 2 ) 5 ' 



V 2 = 



2A(-2A + 3Vl + kr 2 ) 
15(A- Vl + kr 2 ) 2 



(25) 
(26) 



The condition V < 1 is satisfied for A < (18 ~^V l + kr 2 . 

Using Eq. I|24|) and Eqs. we get the following expressions for the energy density, pressure and sound 

velocity of the disk: 



ikA 



A+^l + k{R 2 



P = 



ikA 2 



2tt 
A 



V 2 = 



-A+ ^Jl + k{R 2 +a 2 ) A+ s/i + k(R 2 + a 2 ) 
-A + 2y/l + k(R 2 + a 2 ) 



A - ^/l + k(R 2 + a 2 ) 



(27) 
(28) 

(29) 



The conditions V < 1 and P > are both satisfied if A < + ka 2 . Figure (a)-(d) shows, respectively, a, P, V 
and rotation curves, Eq. (|30[) . as functions of R = R/m for the disk calculated from Buchdahl's solution. 

In figure [7| (a)-(b) we show, respectively, the density p together with pressure p, and sound velocity V of the 
halo along the axis z for A = 0.6; k = 1 and a = 1. Note that in this solution there is no boundary of the fluid sphere: 
the disk is completely immersed in the fluid. 

The tangential velocity v c calculated from metric coefficients l|24|) is 



v? = 



2AkR 2 



I- A/yJl + k{R 2 + a 2 )) {[1 + k(R 2 + a 2 )} 3 / 2 + A[l + k(a 2 - R 2 )}} 



(30) 



For R » a, Eq. (|TT77f> goes as v c = (2A) 1 / 2 /(i? 1 / 2 fc 1 / 4 ). The specific angular momentum follows from Eq. (|2*2*|l and 
Eq. J21: 



V2AkR 2 ( 1 



y/l + k(R 2 +a 2 ) 



[1 + k(R 2 + a 2 )] 1 / 4 



[1 + k(R 2 



-,2\]2 



- AAkR 2 v / l + k(R 2 + a 2 )- A 2 [l + k{a 2 - R 2 )} 



(31) 



B. Narlikar-Patwardhan-Vaidya Solutions la and lb 

Now we shall study the generation of a disk solution with an halo exactly as the one depicted in Fig. [3 We start 
with a solution of the Einstein equations in isotropic coordinates that represents a sphere with radius rt, of perfect fluid 
that on r = will be continuously matched to the vacuum Schwarzschild solution. Narlikar, Patwardhan and Vaidya 
[2~Ij gave the following two exact solutions of the Einstein equations for a static sphere of perfect fluid characterized 
by the metric functions {X,i>i a ) and (A, vw), 
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FIG. 6: (a) The surface energy density a, Eq. 1271 . (b) the pressure P Eq. 128H . (c) the velocity of sound V Eq. 1291 . (d) the 
tangential velocity v c Eq. <!30t for the disk with A — 0.6; k = 1; for a — 1, 2 and 3 as function of R = R/m. 




FIG. 7: (a) The density p and pressure p Eq. Il2ril and (b) the velocity of sound V Eq. Il2(i> for the halo with A = 0.6; k — 
1 and a = 1 along the axis z. 



e x = Cr k , (32) 

e Vl ° = (Ai a r 1 -"+ fc / 2 + B la r 1+n+k / 2 ) 2 for - 2 + v^2 < k < 0, (33) 

e «u _ r V2 + Blb ln ( r )]2 for = _ 2 + \/2, (34) 

where A% a , An, B\ a , Bu, C are constants, and n = \J\ + 2k + k 2 /2. We shall refer to these solutions as NPV la 
and NPV lb, respectively. 
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The density, pressure and sound velocity for for the solutions (A, vi a ) and (A, fib), will be denoted by (p,pi a , Vi a ) 
and (p,Pib, Vib), respectively. We find, 



P 

Pla 
Plb 

v? 

v la 



-k(k + A)r 



-2-k 



32ttC 

= 327rCr2+k( l a + Biar2n) {Ma + 8(1 - n) - 4fc(n - 3)] 

-B la [ik 2 + 8(1 + n) + 4fc(n + 3)] r 2 ™} , 

Aib + B lb ln(r) + gyggib 
16nC[A lb + flibln(r)] ' 

= fc(fc + 4)( ^ flinr2m)2 [-3^ + 8(n - 1) + 4 fc (n - 3)] 

-B 2 a r An [3k 2 + 8(1 + n) + Ak{n + 3)] - 2A la B la r 2n [3fc(fc + 4) + 8(1 - n 2 )] } , 
2B\ h r^ 2 



[A lb + B lb ]n(r)]< 



(35) 

(36) 
(37) 

(38) 
(39) 



The condition of continuity of the metric functions (A, v) given by (|32|I - H34(I and the corresponding functions 
in Eq. (0 at the boundary r — r b leads to following expressions: 



m 
2n 



C = 



Ala = 



3fc 2 + 8(l + n)+4fc(n + 3) 



16nr 



l-n+fe/2 



Bla = 



k + A / 

3fc 2 + 8(l-7i)-4/c(n-3) 



16m' 



l+n+fc/2 



2V2 + ln(r b ) 
Alb = TT7^ ) — 



1 



Ar 



1/V2 



Ar 



1/V2' 



(40) 
(41) 

(42) 



Vi a has its maximum at r = 0, and Vib at r = r b . Condition Vi b (r b ) < 1 is satisfied if r b < A 1 /^ 2 . 

Using Eq. H32(l - (|34|l in Eq. we get the following expressions for the energy density, pressure and sound 

velocity of the disk: 



ka 



Pla = 



AttVCK^/ 4 ' 

a [A la (k - n + 1) + Bia(k + n + l)H n ] 



Plb = 



[Ala + Bi a K n ] 
2Ai b 



2-1) 



2-Bit 



V? 

v la 



V 2 
v lb 



4 7 rVC^ 1 /2+v^/4[ 2 ^ lb + Bib ln(ft)] 
+Bib(V2 - 1) ln(ft) 

= - 7 \-Bf a TZ n (k 2 + 5k + nk + An + 4) 

k(k + 4) [Ai a n- n / 2 + Bi a n n / 2 } 

+A 2 la K- n {-k 2 - 5k + nk + An - 4) + 2Ai aJ Bi a (-/c 2 - 5k + An 2 - A)] , 

= 9^9 a 1 p i 7^2 k 2 fc \/2 In 2 (ft) + 2Bib(2B u + B U V2 + 2A U V2 
2(2Aib + Bibin(TZ)y L 

x ln(ft) + 4A lb B lo (2 + V2)+ AV2A 2 b + 8B 2 



(43) 
(44) 

(45) 

(46) 
(47) 



where 1Z — R 2 + a 2 . V\ a and Vi b have their maximum values at R = 0. Because the expressions are rather involved, 
the restrictions on the constants to ensure the velocities are positive and less then one is best done graphically. The 
curves of er, P and V as function of R = R/m with parameters k = —1/2; r b = 2 and a = 0.5, 1.0,1.5 are displayed 
in figure [S] (a)-(c), respectively. Figure ED (a)-(b) show the density p, pressure p and velocity of sound V for the halo 
with parameters k = —1/2; r b — 2, for a — 0.5 along the axis z. The same physical quantities are shown in figures 
^2 ( a )~( c ) and^|(a)-(b) with k = —2 + \/2. We note that a and P are continuous at the boundary between the 
internal and external parts of the disk, but the velocity of sound has a discontinuity. 
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FIG. 9: (a) The density p Eq. I|85|l and pressure p Eq. . (b) the velocity of sound V Eq. I|38|l for the halo with k = —1/2; r;, = 
2; for a = 0.5 along the axis z. 



The tangential velocity v c is given by 



2 A ln (l - n + fc/2) + gOl + n + fc/2)QB 2 + a 2 f 
" ! " [A la + i? la (i? 2 + a 2 )"][i? 2 (fc + 2) + 2a 2 ] 

, ., 2\/2A lfc + 5 lb [4 + V21n(i? 2 + a 2 )] 



clb [2A lfc + 5 lb ln(i? 2 +a 2 )][2 a 2 + V2i? 2 ]' 
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FIG. 10: (a) The tangential velocity v c i Eq. igBJ and (b) the curves of hj~ Eq. JJJJ witn k = -1/2; n = 2; for a = 
0.5, 1.0, and 1.5 as function of R = R/m. A region of instability appears on the disk generated with parameter a = 0.5. 
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FIG. 11: (a) The surface energy density a Eq. 1430 . (b) the pressure P Eq. I|45|l . (c) the velocity of sound V Eq. I|470 for the 
disk with k = — 2 + \/2; rj, = 2; for a — 0.5, 1.0 and 1.5 as function of 7? = R/m. 



and the specific angular momentum h, 

h la = VCR 2 (R 2 + a 2 ) k /\ 



! A la (l -n + k/2) + B la (l + n + k/2)(R 2 + a 2 )" 
A la {a 2 + nR 2 ) + B la {a 2 - nR 2 )(R 2 + a 2 ) n ' 



V 2a 2 [A lb + S lb ln(Vi? 2 + a 2 )] - 2B lb R 2 



(50) 



(51) 



In figure ITUI (a)-(b), the curves of tangential velocity Eq. H48fl an d h~ Eq. 1)51) |). respectively, are displayed as 
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FIG. 12: (a) The density p Eq. and pressure p Eq. I|37|l . (b) the velocity of sound V Eq. Ij39^ for the halo with k 
—2 + \/2; r b — 2; for a = 0.5 along the axis z. 
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FIG. 13: (a) The tangential velocity v c2 Eq. (gHJ, 0) the curves of hj~ Eq. f3TJ with fc = -2 + ^/2; r b = 2; for a 
0.5, 1.0, and 1.5 as function of R = R/m. As in the previous case, the same region of instability occurs. 



functions of R = R/m with k = —1/2; rj, = 2; a = 0.5, 1.0,1.5. The same quantities are shown in figure IT^T a'l-fb) 
with k = —2 + For a = 0.5 the disks have a small region of unstable orbits immediately after the "boundary 
radius" . 



C. Narlikar-Patwardhan-Vaidya Solutions 2a and 2b 



Like in the previous subsections we study the generation of a disk solution with an halo exactly as the one 
depicted in Fig. |2 We also start with a solution of the Einstein equations in isotropic coordinates that represents a 
sphere of radius rj, of perfect fluid that on r — r\, will the continuously matched to the vacuum Schwarzschild solution. 
We will use another two solutions found by Narlikar, Patwardhan and Vaidya [24| that we shall refer as NPV 2a and 
NPV 2b, respectivelly. that are characterized by the metric functions (A, Via) and (A, ^2h), 



(Air 1 +"/ 2 +A 2 r 1 -"/ 2 ) 2 

( Air i+ n /2 +A2r i- n/ 2y 

_ [B lb +B 2b \n(r)] 2 

~ (A ir i/^ + i !r -i/V2) 2 



for V2 < n < 2, 
for n = V2. 



(52) 

(53) 
(54) 



where the A's and B's are constants and x = \ / 2n 2 — 4. The solution (A, V2 a ) with n = 2 corresponds to Schwarzschild's 
internal solution in isotropic coordinates (see for instance Ref. |22jV This solution has constant density and is 
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conformally flat when B\ a = 0. 

The density, pressure and sound velocity for the solutions (A, V2a) & n d (A, V2b), will be denoted by (p,P2a, ^2a) 
and {p,P2b,V2b), respectively. We find, 



P = 

P2a 

■ 

P2b = 



(4 - n 2 ) (Air™/ 2 + A 2 r-"/ 2 ) 2 + l2n 2 A x A 2 
-12n 2 A 1 A 2 + (3n 2 - 4) (A x r n/2 + A 2 r- n/2 ^j 



1 

32^ 
1 

32^ 

2nx(B 2a - B la r x ) 



B 2a + B la r* 
1 

167r(Bi 6 + B 2b Hr)) 



(A 2 r n - A\r- n ) 



(B lb + B 2b Hr)) (Alr^ + Alr 



10AM; 



v. 



+2V2B 2b [Alr-^ 2 - A\r 



1 



2„,V2 



{B\ a [A 2 r"(3n 2 + 2nx - 4) 



V. 



~ (4 - n 2 ) (Ajr n - A 2 r - n ) {B 2a + B la r x ) 2 
+A 2 2 r- n {-3n 2 + 2nx + 4)] + B 2 la r 2x [A 2 r"(3n 2 - 2nx - 4) 
+A 2 r-"(-3n 2 - 2nx + 4)] + 2B la B 2a r x [A 2 r"(3n 2 - 2x 2 - 4) 
+A 2 r~ n (-3n 2 + 2 ; r 2 + 4)]}, 

' [a\t^ 2 \{B lb + B 2b Hr)) 



2b 



[B lb + B 2b ln(r)] 2 (a 2 t^ - A 2 r -^ 



-A 



2 -V2 



: (B lb + B 2b ln(r) - 2V2B 2b ) + 2B- 
(B lb + B 2b ln(r))(B lfc + B 2b ln(r) + 2V2B 2b ) + 2£ 2 2 b ] } 



(55) 



(56) 



(57) 



(58) 



(59) 



The condition of continuity of the metric functions (A, v) given by i|52 [l — 1|54| ) and the corresponding functions 
in Eq. © at the boundary r = r b leads to following expressions: 



A x 



An 



Si 



So 



B 



lb 



B 



2b 



2+n/2 /-. 
nT b \ l ' 2r b J 



m In mn 

fb\ 1 

2 V 2 4r fc 



2-n/2 



2/7 



'•ft 



-4r 2 



4r 2 



(l - 2s) - to 2 + 2xr 2 (l - ^ 

4-r /2 (l + ^) 4 
(l-^)+^ 2 + 2.^ 2 (l~f0 

^5~ S/a (! + ^) 4 



11m 

2 n 4r b 



^ (l + ft) 
77i 2 — 8mrb + 4r 2 



4r 2 - to 2 + (to 2 - 8TOr b + 4r 2 ) ln(r fc )] 



(60) 
(61) 

(62) 

(63) 
(64) 
(65) 



Y 2a has its maximum at r = r b , and V26 at r = 0. 

Using Eq. (|52|I - H54() in Eq. 1(311-©! we get the expressions for the energy density, pressure and sound velocity 
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of the disk: 



4tt 



At(2 + n)K- 1/2+n/i + A 2 {2 - n)K- 1/2 - n/4 
a 



(66) 



Pla = 7 ttz — r. ; 7-t B la A 1 (2 + 2n - xYR. 

+B la A 2 (2 - 2n- a;)^ (;r - n)/4 + B 2o Ai(2 + 2n + x )T^- x+n ^ 
+5 2QJ 4 2 (2-2n + a ;)ft- (a: + rl H 



(m+n)/4 



(67) 



P 2h = - 



{2(1 + V2)B lb AiK- 1/i+VS/ * 

B 2b A 1 n- 1 ' 2+ ^/ i 



V 2 

> 9„ 



47r[25i 6 + 5 2b ln(ft)] 
2(1 - V2)B lfc A 2 ^- 1/2 -^ /4 + [(1 + \/2) ln(ft) - 2 

(1 - V2) fn(7£) - 2 

1 



i^A^- 1 / 2 -^/ 4 } 



(68) 



/2 



2(4 - n 2 ) [A 2 + A 1 U n / 2 ] [B 2a + B la K*/*]' 
[Bf a (n-2){2n-x + 2)K x +B 2 a {n-2)(2n + x + 2) 

-4B la B 2a (-n 2 + n + x 2 + 2)TZ x/2 ] + A 2 [B 2 la (n + 2){2n + x- 2)TZ X 
+B 2 a (n + 2)(2n - x - 2) - 4B la B 2a (~n 2 -n + x 2 + 2)7T /2 ] } , 



V& = - 



2B lb + 2B 2b + B 2b ln(K) 



{ V21n(^)-4 



2[2Bi b + B 2b \n(TZ)] 2 [A 2 + AxTZ 1 '^] 
B^AxV}!^ - [V21ii(72.) + 4] B 2b A 2 + 2V2B lb ^TZ 1 ^ - A 2 ] } 



(69) 



(70) 



where Tl = R 2 + a 2 . The curves of a, P and V as function of R = R/m with parameters n — 1.8; m = 0.5; r b = 
2 and a = 0.5, 1.0,1.5 are displayed in figure ITTI (a) - (c), respectively. Figure ITSI (a) - (b) shows the density p, 
pressure p and velocity of sound V for the halo with parameters n = 1.8; m = 0.5; r b = 2, for a — 0.5 along the axis 
z. The same physical quantities are shown in figures [T7|and EH with n = V2. 
The tangential velocity v c is given by 



"c2o 



R 



2 A 2 [B la (n + x)K x ! 2 + B 2a (n - x)} - A^^jB^n - x)K x ' 2 + B 2a (n + x)] 
[B la K x / 2 + B 2a ][A 1 TZ n / 2 (2a 2 - nR 2 ) + A 2 (2a 2 + nR 2 )} 



2 _ D2 A 2 [2V2B lb + B 2b (4 + V2\n{1l))} - A^^V^Bu + B 2b (-4 + V2ln(1l))} 

S c 2b — K 



[2B lb + B 2b \n{'R)][A 1 'R^ 2 / 2 {2a 2 - ^/2R 2 ) + A 2 (2a 2 + V2R 2 )} 
and the specific angular momentum h 

R 2 n -l/2+n/4 



(71) 
(72) 



l'2a 



h 2b = 



{AiTZ n / 2 + A 2 fl 2 ^B la (2a 2 - xR 2 )K x l 2 + B 2a (2a 2 + xR 2 ) 

[A 2 [B la {n + x)n x ' 2 + B 2a (n - x)] - A{R n ' 2 [B la (n - x)Tl x ' 2 + B 2a (n + x)]j 

^-1/2+72/4 

[AJZ^/ 2 + A 2 ]V 2 y/4tfB lb + 2B 2b [-2R 2 + a 2 ln(ft)] X 
[A 2 [2V2B lb + B 2b {4 + y/2hx{H))] - A 1 1l" /2/2 [2V2B lb + B 2b {-4 + V2hx{K))]\ 



1/2 



1/2 



(73) 



(74) 



In figure [TBI (a)-(b), the curves of tangential velocities Eq. H71|) and Eq. I|73|l . respectively, are displayed as 
functions of R — R/m with n = 1.8; m = 0.5; r b = 2; a = 0.5, 1.0,1.5. Figure ITSI shows the same quantities with 
n = \2. Unlike solution 2, no unstable circular orbits are present for the disks constructed with these parameters. 
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FIG. 14: (a) The surface energy density a Eq. Ititjl . (b) the pressure P Eq. 16711 . (c) the velocity of sound V Eq. 1691 for the 
disk with n = 1.8; m = 0.5; rt = 2; for a — 0.5, 1.0 and 1.5 as function of R = R/m. 
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FIG. 15: (a) The density p Eq. 1551 and pressure p Eq. I56H . (b) the velocity of sound V Eq. 1581 for the halo with n = 1.8; m 
0.5; rb = 2; for a — 0.5 along the axis z. 



V. DISKS WITH COMPOSITE HALOS FROM SPHERICAL SOLUTIONS 



We study two examples of disks with halos constructed from spheres of fluids with two layers as the ones 
depicted in Fig. 
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a = 0.5 
a- 1.0 
a- 1.5 



FIG. 16: (a) The tangential velocity v ca Eq. ED, ( b ) tn ^ curves of hf± Eq. EHJ with n = 1-8; m = 0.5; r b = 2; for a 
0.5, 1.0, and 1.5 as function of R = R/m. The disks have no unstable orbits for these parameters. 
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FIG. 17: (a) The surface energy density a Eq. 1)660 . (b) the pressure P Eq. I|fi8|l . (c ) the velocity of sound V Eq. I|7O0 for the 
disk with n — y/2; rt = 2; m — 0.5; for a — 0.5, 1.0 and 1.5 as function of R = R/m. 



A. Internal Schwarzschild solution and Buchdahl solution 



Let us consider a fluid sphere is formed by two layers: The internal layer, < r < ri, will be taken as the 
internal Schwarzschild solution (solution 2a with n=2), 

*-W+b,)' eA = 1 

(A^ + A,)' {A,r' + A 2 f 
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FIG. 18: (a) The density p Eq. II55|I and pressure p Eq. \F>7}. (b) the velocity of sound V Eq. 1)590 for the halo with n = \/2; rt = 
2; for a = 0.5 along the axis z. 
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FIG. 19: (a) The tangential velocity v cb Eq. JZU, 0) The curves of Eq. CD with n = ^2; m = 0.5; r 6 = 2; for a 
0.5, 1.0, and 1.5 as function of R = R/m. The disks have no unstable orbits for these parameters. 



The external layer, r > r\ is taken as the Buchdahl solution, 



1 - 



c 



\/l+kr 2 



1 



C 



1 



c 



VT + kr 2 



Note that the external layer has no boundary, i.e., this layer has infinite radius 



(76) 
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FIG. 20: (a) The surface energy density a, (b) the pressure P, (c ) the velocity of sound V for the disk generated from spherical 
fluid layers Eq. H751 - 176H with m = 1; k = 1; ri = 2; for a = 0.5, 1.0 and 1.5 as function of 7? = R/m. 



According to the continuity conditions at r = ri, the constants are related through: 



A 1 



D 



Bo 



Ck 



(c+y/l + kr\ 
Ck 



(1 + krf) 1 



A, = 



4- g. 

T (l+fcrf) 3 / 2 



C(l - 3fcr2) - ^ 2 (1 - fcri) + 2(1 + kr\f/ 2 



1 - 



(1 + kr\Y + 2C v /l + fc?'2 + C 2 (l - kr\) 
c 

Bui. 



(77) 
(78) 

(79) 
(80) 



With these relations, one verifies that, using Eq. (|2"7|) and Eq. ljrjfj|) . Eq. (f^5|) and Eq. ljrj7|) . both the energy density 
and the pressure are continuous at the radius R = y/ r\ — a 2 of the disk. 

Figure [3U| (a)-(c) shows, respectively, c, P and V for the disk obtained from fluid layers Eq. (|75() - (|76|l with 
parameters m = 1; k = 1; r\ = 2; for a — 0.5, 1.0 and 1.5 as function of R = R/m. The density p 1 pressure p and 
velocity of sound V for the halo along the z axis with the same parameters for a = 0.5 is shown in figure PHI 

In figure (a)-(b), the curves of tangential velocities and of h^j~, respectively, are displayed as functions of 

R = R/m. 
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FIG. 21: (a) The density p and pressure p, (b) the velocity of sound V for the halo formed by fluid layers Eq. I|75p - I|7ti[l with 



; n = 2 and a 


= 0.5 along the axis z. 
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FIG. 22: (a) The tangential velocity v c , (b) the curves of h~ for the disk generated from fluid layers Eq. I751 - I76H with 

m = 1; k = 1; n = 2; for a = 0.5, 1.0 and 1.5 as function of R — R/m. The disks obtained with these parameters have no 
unstable orbits. 



B. NPV Solution 2b with n = \[2 and NPV solution lb with k = -2 + \[2 

Now we consider a sphere composed with two finite layers: The internal layer, < r < r 1; is taken as the NPV 
solution 2b with n — \/2, 



(Si + £ 2 ln(r)) s 



1 



(81) 



The external layer, n < r < r%, is taken as the NPV solution lb with k = —2 + \/2, 

e v = (A 3 + B 3 ln(r)) 2 , e x = Cr- 2+V *. (82) 
The spacetime outside the sphere, r > r 2 , will be taken as the Schwarzschild's vacuum solution is isotropic coordinates. 



(1 + f 



2 ■ 



1+ 27 



(83) 
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a=0.9 



FIG. 23: (a) The surface energy density a, (b) the pressure P, (c ) the velocity of sound V for the disk generated from spherical 
fluid layers Eq. 1811 — 1821 with n = 1; r<z = 2; for a — 0.3, 0.6 and 0.9 as function of 7? = R/m. 



In this case the pressure should be zero at r = r%. The continuity conditions at r = r\ and r = r-i give the relations: 



m 

r-2 



V2(2 - V2) 
l + \/2 



64r. 



2-V2 



(1 + V2) 4 



, B 3 



4r, 



A, 



2V2 + ln(r 2 ) 



Ai = 0, A 2 = -=, B 2 = -^= 

Vc 2V2C 



,%/2 



2V2C 



(A 3 + B 3 ln(n)) (1 - r^) + 2V2B E 



(A 3 + B 3 ln(n)) (2V2r-^ - Info) + - 2\7iB 3 ln(n) 



(84) 

(85) 
(86) 



Using Eq. and Eq. (JSSJ), the energy density of the disk at R = y/r'f — a 2 is continuous, but not the pressure. The 
difference between Eq. (|68|) and Eq. H45J1 is 



AP 



a{r\ 



v / 2/2 -V2/2n 



167rVCri(A 3 + B 3 \n(n)) 

„\/2/2 -n/I/2 



V2(A 3 + B 3 ln(n)) + 4B 3 



(87) 



The pressure is continuous if r?*'* — r, v ^ 7 ^ = — > ri = 1. 

Figure [2S1 ( a )~( c ) shows, respectively, a, P and U for the disk obtained from fluid layers (|81|I - H82() with 
parameters r\ = 1; r-i = 2; for a — 0.3, 0.6 and 0.9 as function of R = R/m. The density p, pressure p and velocity 
of sound V for the halo along the z axis with the same parameters for a — 0.3 is shown in figure l2*4l 

In figure 1251 (a)-(b), the curves of tangential velocities and of h~, respectively are displayed as functions of 

R = R/m. In this case, regions of unstable orbits exist for parameters a = 0.3 and a = 0.6. 
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FIG. 24: (a) The density p and pressure p, (b) the velocity of sound V for the halo formed by fluid layers Eq. I|81^ - I|82|l with 
Ti = 1; ra = 2 and a = 0.3 along the axis z. 




FIG. 25: (a) The tangential velocity v c , (b) the curves of h~ for the disk generated from spherical fluid layers Eq. H8H 18211 
with n = 1; V2 = 2; for a — 0.3, 0.6 and 0.9 as function of R = R/m. Regions of unstable circular orbits appear for the disks 
obtained with parameters a — 0.3 and 0.6. 



VI. DISCUSSION 



The "displace, cut and reflect" method applied to solutions of Einstein field equations in isotropic coordinates 
can generate disks with positive energy density and equal radial and azimuthal pressures (perfect fluid) . With solutions 
of static spheres of perfect fluid it is possible to construct disks of perfect fluid surrounded also by perfect fluid matter. 
As far we know these are the first disk models of this kind in the literature. 

All disks constructed as examples have some common features: surface energy density and pressures decrease 
monotonically and rapidly with radius. As the "cut" parameter a decreases, the disks become more relativistic, with 
surface energy density and pressure more concentrated near the center. Also regions of unstable circular orbits are 
more likely to appear for high relativistic disks. Parameters can be chosen so that the sound velocity in the fluid and 
the tangential velocity of test particles in circular motion are less then the velocity of light. This tangential velocity 
first increases with radius and reaches a maximum. Then, for large radii, it decreases as 1/y/R, in case of disks 
generated from Schwarzschild and Buchdahl's solutions. The sound velocity is also a decreasing function of radius, 
except in solution NPV 2a with \[2 < n < 2, where it reaches its maximum value at the boundary. In principle other 
solutions of static spheres of perfect fluid could be used to generate other disk + halo configurations, but it is not 
guaranteed that the disks will have the characteristics of normal fluid matter. 

We believe that the presented disks can be used to describe more realistic model of galaxies than most of the 
already studied disks since the counter rotation hypothesis is not needed to have a stable configuration. 

We want to finish our discussion by presenting a table that summarizes our results about disks in an unify 
manner. 
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Disks properties 



Solution 


metric 


matching 


energy 


pressure 


sound 


angular 




coefficients 


conditions 


density 




velocity 


momentum 


~r~y , l /-^ i I'll 

External bcnwarzscnud 


© 




cm 


CD 


C3 




Buchdahl 






<E3 




(E3 


<E3 


NPV la 


02l G3J 




63 


62) 


63 




NPV lb 


(EH, p3| 


1001), EH 


63 


63 


63 


l r 1 1 

<E3 


NPV 2a 


ij52l-Oi 




(HSJi 


E3 


(E3 


<E3 


NPV 2b 








(E3 


CZ3 


G3 


Internal Schwarzschild 














+Buchdahl 




CZ3-JE3 




123 


JMJ, (E3 


CZ3, EH) 


NPV 2b+NPV lb 


IHU, JB2J), 


(JElJi-iJHii 


JSU, 63, 


(E3, 63, 


dZ3, 63, 


(13, EH), 


+external Schwarzschild 


(ESI 




GDI 


CD 


C3 





In the table we list the seed metric coefficients, matching conditions at the boundaries and relevant physical 
quantities of all disks studied in this work. The numbers refer to the equations presented along the paper and NPV 
stands for Narlikar, Patwardhan and Vaidya as before. 
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